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ABSTRACT 

Models of electroweak symmetry breaking with more than a single doublet 
of Higgs scalars contain a neutral pseudoscalar boson. The production of such 
a pseudoscalar in hadron collisions proceeds primarily via gluon fusion through 
a top-quark loop (except for those models in which the pseudoscalar coupling to 
bottom quarks is strongly enhanced). We compute the QCD corrections to this 
process in the heavy-quark limit, using an effective Lagrangian derived from the 
axial anomaly. 
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The primary mission of the next generation of colhders is the study of the 
physics of electroweak symmetry breaking. Since the Higgs mechanism provides the 
most promising scenario for the symmetry breaking, the search for the particles of 
the Higgs sector is of primary importance. The minimal model of electroweak sym- 
metry breaking contains one complex scalar doublet, three components of which 
become the longitudinal degrees of freedom of the and Z. The remaining com- 
ponent of the doublet is the so-called Higgs boson. However, the minimal model 
with one doublet has no a priori justification and there are several motivations for 
considering models with enlarged Higgs sectors, either containing more doublets, 
singlets or more exotic representations. For example, supersymmetric models re- 
quire at least two doublets. Likewise, at least two doublets are required to produce 
CP violation in the Higgs sector[l]. 

Models with enlarged Higgs sectors have a richer particle content than the 
minimal model; in general, neutral pseudoscalars (with respect to their fermion 
couplings) and charged scalars as well as extra neutral scalars are present. In this 
letter we will study the QCD corrections to the production of a Higgs pseudoscalar 
{A'^) via gluon fusion. This process proceeds primarily through a top-quark loop 
(unless the coupling to bottom quarks is greatly enhanced). We will focus on the 
case of a light pseudoscalar and work in the heavy- top- quark limit : mt ^ Ma. 

To fix our normalization we take the coupling of the pseudoscalar to top quarks 

to be mt'j5/v where v = 246 GeV[2]. The lowest-order amplitude for gg A'^ is 
well known and takes the form [3] 

M{g{kl),g{kl) ^ A') = -igAS'''e^'''''hpk2aTf{T), (1) 
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with r = Ami /Ml, 9 A = OLsjl'nv, and 



sin-i ( Jljr 



if r > 1, 



-1/4 [ln(r;+/r;_) - Z7r]2 , if r < 1, 



where r\± — (1 ± \J\ — r). We see that in the large mt hmit the amphtude for 
gg — > is independent of the top quark mass, just as in the scalar case: T/(r) — > 1 
as T — > oo. The heavy-top approximation of the amplitude is accurate to within 
5% for > 2Mf and to within 10% for mf > M^. 

The amplitude in Eq. (1) can be computed by evaluating the triangle diagram 
with a top quark in the loop and taking the limit — > oo or instead by noticing 
that this amplitude is related to the axial anomaly[4]. Let js = ip'ysip be the axial 
current and = tpij^jsip be the axial vector current. The anomaly equation 
reads: 

d,j^ = -2mtj, + (2) 

where G^j^ is the field-strength tensor for SU(3) and G^j^ is its dual, G^j^ = 
eMi'po-(7o^_ In the heavy quark limit the left side of Eq. (2) vanishes; it is pro- 
portional to more powers of external momenta than the other terms [5]. Therefore, 
in the heavy quark limit the matrix element of between gluon states is given 

2mt{g\Mg) = -^{9K,G;jg). (3) 

Equivalently, one may treat the interactions of gluons with the A'^ in the heavy- 
quark limit as arising from the effective Lagrangian[6] 



^eff = ^G«^G«^>1°. (4) 



The power of Eq. (4) for our problem comes from the Adler-Bardeen theorem [7] 



which states that Eq. (2) is true to all orders in perturbation theory. Therefore, Eq. 
(4) provides the correct effective Lagrangian with which to compute radiative cor- 
rections. The advantage is that amplitudes which correspond to two-loop diagrams 
in the original theory are one-loop diagrams in the effective theory. A similar La- 
grangian can be written for the interaction of photons with the pseudoscalar[8,9]; 
in that case one sees immediately that the 0{as) corrections to 77 vanish 
in the mt ^ 00 limit. 

The effective Lagrangian, Eq. (4), leads to two- and three-gluon vertices with 
the (as well as a four-gluon vertex which is irrelevant to the current study): 



where the ki are the gluon momenta directed inward. 

We use dimensional regularization in computing the radiative corrections. How- 
ever, the e-tensors in Eq. (5) are intrinsically four-dimensional objects and must 
be treated as such. The product of two e-tensors can be written in terms of g'^''' 
the metric tensors in the four-dimensional sub-space: 



where g^^gjiv — 4. For simplicity, we take the incoming momenta to be in four 
dimensions (this is simply a choice of frame). The lowest order cross section is 
then, averaged over colors and polarizations in n = 4 — 2e dimensions. 



V4(kf., kl, kL) = ggAr"'e^'''"'(h + k2 + ks) 



(5) 
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where z — M\l s and s — {k\^ kif'. 
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The virtual radiative corrections are given by the diagrams in Figure 1, where 
the amphtude involving the four-gluon vertex vanishes due to the antisymmetry of 
the ggA^-vertex. The other diagram yields the result for the virtual corrections 



where 

'47r/i2 



77 = r(l + e) 



M2 



A 



The real corrections are given by the diagrams in Figure 2. The amplitudes 
squared for the various processes are, averaged over colors and spins, 



\M{gg - = 87V,«,a(,) I — ^ - 2es ( 1 + 2^3 • hs' 



{ 2 2 2 

-^^+2^^3-^3|, 



(9) 

where s, and u are the familiar Mandelstam variables and is the momentum 
of the final state quark or gluon, with the hat denoting the (n — 4)-dimensional 
components. Since the initial state particles are taken to be in four dimensions 
/c3 • /c3 is the only quantity which depends on the {n — 4) -dimensional components. 
Under integration over the angular variables transverse to the incoming particles 
we find [10] 

j dn-rh ■ h = -|^3^^T j dn^, (10) 

where fcx is the transverse momentum of the outgoing gluon. So the terms with 
ks ■ ks will contribute terms of C(l), but only when multiplied by double poles, 
or u~'^ (this was used in simplifying Eq. (9)). 



Integration over n-dimensional phase space gives the cross sections 



a{gg ^ gA^) 



OisNcCr^e) 



TT 



1 



TT 



2z 



z 



+ 2 



(1 - ^)+ 

\4 , 4 



H h-2(l - z) 

z 



'^{Q9 - = ^^^(1 - (-^ + 21og(l - [1 + (1 - 



--(l-z){7-3z)}, 



where 



2asC^ao 



1 - zf, 



(11) 



It is interesting to note that the real corrections, when written in terms of the 
lowest-order cross section are identical in form to those for the scalar case [6]. 

We see that the terms of (9(l/e^) cancel between the real and virtual diagrams. 
The terms of (9(l/e) may be absorbed into redefinitions of the parton distribution 
functions in the usual factorization procedure. We use the MS prescription which 
fixes the subtraction terms as 



a 



4p = ^%)2^P,,(^)(47r)T(l + e)-, 
4p = ^^(e)^^P5(^)(4vr)T(l + e)l 



(12) 



where the splitting functions are defined 

z 1 — z 



Pgg(z) = 2Ar, 



l-z)+ 
[1 + (1 - z)'] 



+ z{l - z) 



+ PoS{l - z) 



(13) 



where (3q — ^-Nc — ^Nf with Nf being the number of fight quark fiavors. The 



charge renormalization subtraction term is 

cxeh = -^(47r)T(l + e)a(,)<5(l - z)^. (14) 



When all the contributions are included the final results are 




+ z{l-z) 



(15) 

Since the real corrections in our result are the same as for the scalar case, the 
only difference between the scalar and pseudoscalar results is the coefficient of the 
(5(1 — z) term coming from the virtual result: Nc{7r'^/3 + 2) for the pseudoscalar 
and {Ncn'^/3 + 5/2Nc — 3/2Cp) for the scalar. The numerical similarity of the two 
constant terms [an accident of SU(3)] means that the ratio between the next-to- 
leading order result and the lowest-order result, or the 'K- factor,' will be almost 
equal for the two processes. 

The radiatively corrected cross section is plotted for proton-proton collisions 
at VS — 40 TeV and = 15.4 TeV in Figure 3 and Figure 4 for several choices 
of renormalization scales and using HMRSB parton distributions [11]. The ratio 
between the radiatively corrected result and the lowest order result (computed 
with the two-loop ag) ranges from about 2.6 for Ma = 50 GeV to about 2.2 for 
Ma — 200 GeV for both values of v^- As is the case for the scalar Higgs, the 
contribution from the gg initial state dominates with about half the correction 
coming from the 5(1 — z) term. 
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FIGURE CAPTIONS 

1) Feynman diagrams for the virtual corrections to gg ^ A^. 

2) Feynman diagrams for the processes a) gg — > gA^ and b) qq — > gA^. The 
diagram for qg — > qA'^ is a crossing of diagram b. 

3) Radiatively corrected results (utot) for A*-^ production in proton-proton col- 
hsions at the SSC, y/S — 40 TeV, with parton distributions from Ref. 11 
for several values of the factorization/renormalization scale fi. Shown for 
comparison is the lowest order cross section (computed with the two-loop 
as) with II — Mj^. 

4) Same as Figure 3 for the LHC, y/S = 15.4 TeV. 
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